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Abstract 

In the simulation of intermediate energy heavy ion collisions by using the isospin 
dependent quantum molecular dynamics, the isospin effect on the process of multi- 
fragmentation and dissipation has been studied. It is found that the multiplicity of 
intermediate mass fragments Ni m f for the neutron-poor colliding system is always 
larger than that for the neutron-rich system, while the quadrupole of single par- 
ticle momentum distribution Q zz for the neutron-poor colliding system is smaller 
than that of the neutron-rich system for all projectile-target combinations studied 
at the beam energies from about 50MeV/nucleon to 150MeV/nucleon. Since Q zz 
depends strongly on isospin dependence of in-medium nucleon-nucleon cross section 
and weakly on symmetry potential at the above beam energies, it may serve as 
a good probe to extract the information on the in-medium nucleon-nucleon cross 
section. The correlation between the multiplicity Ni m j of intermediate mass frag- 
ments and the total numer of charged particles iV c has the behavior similar to Q zz , 
which can be used as a complementary probe to the in-medium nucleon-nucleon 
cross section. 
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1 Introduction 



In recent years, the establishment of radioactive beam facilities at many labora- 
tories over the world and the use of radioactive beams with large neutron or proton 
excess have offered an excellent opportunity to investigate the isospin-dependence of 
heavy ion collision (HIC) dynamics^ 1-4 ]. This kind of study has made it possible to ob- 
tain the information about the equation of state (EOS) of asymmetric nuclear matter 
ranging from symmetric nuclear matter to pure neutron matter and the information 
on isospin-dependence of in-medium nucleon-nucleon (N-N) cross section, which are 
significantly important not only in understanding nuclear properties but also in ex- 
ploring the explosion mechanism of supernova and the cooling rate of neutron stars. 
However two essential ingredients in HIC dynamics, the symmetry potential of the 
mean field and the isospin-dependent in-medium N-N cross section, have not been 
well determined so far. Recently, Li et al.' 1,5 ' 6 ' made use of the isospin dependent 
transport theory to investigate nuclear symmetry energy and showed that the rate of 
pre-equilibrium neutron-proton emitted in intermediate energy HIC is sensitive to the 
density dependence of nuclear symmetry potential, but insensitive to the incompress- 
ibility of symmetric nuclear matter and the in-medium N-N cross section. Also, R.Pak 
and Bao-An Li et al. have found that the isospin dependence of collective flow and 
balance energy are mainly originated from the isospin-dependent in-medium N-N cross 
section' 5 ^ 7-11 ]. Recently we have found that nuclear stopping can be used as a good 
probe for exploring the in-medium N-N cross section in HIC in the beam energy rang- 
ing from above Fermi energy to 150MeV/nucleon. However, it is still not clear how 
the nuclear stopping depends on the neutron-proton ratio of the colliding system^ 12 !. 
In viewing that little information is known about the in-medium N-N cross section and 
its isospin dependence up to now, it is thus very desirable to find an efficient way to 
gain such kind of knowledge. 

It is known that many intermediate mass fragments(IMF) are emitted in the process 
of intermediate energy HIC and that the element distribution and fragment multiplicity 
indicate a strong correlation between the multiplicity of intermediate mass fragments 
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Ni m f and the total number of charged particles N c The question is whether the 
multifragmentation and dissipation, especially the N im f and Q zz are sensitive to the 
neutron-proton ratio of the colliding system. To answer this question, we have investi- 
gated the isospin effect on the process of multifragmentation and dissipation in HIC for 
the colliding systems with different ratios of neutron to proton by using isospin depen- 
dent quantum molecular dynamics (IQMD). To increase the efficiency of detectors and 
the statistics of N im f and Q zz , the reversible kinetic reactions with heavy projectiles on 
light targets are suggested to have more intermediate mass fragments emitted towards 
forward angles. 

The calculated results show prominent isospin effect for the multifragmentation 
Ni m f and the quadrupole of single particle momentum distribution Q zz for the collid- 
ing systems with different neutron or proton excesses. The multiplicity of intermediate 
mass fragments N im f of neutron-poor systems is always larger than that of neutron- 
rich systems, while the quadrupole of single particle momentum distribution Q zz for 
neutron-poor systems is smaller than that of neutron-rich systems at the beam ener- 
gies from 50MeV/nucleon to f50MeV/nucleon for all the reversible colliding systems 
studied here. The mechanism responsible for the above isospin effect can be found 
from the fact that the mean N-N cross section for a neutron-poor system is larger than 
that of the corresponding neutron-rich system with the same masses of projectile and 
target. The calculated results also show that the correlation between the multiplicity 
of intermediate mass fragments Ni m j and the total number of charged particles N c , de- 
pends strongly on the isospin dependence of in-medium N-N cross section and weakly 
on the symmetry potential in the chosen beam energy region. In Ref. [12], it is found 
that Q zz is very sensitive to the isospin dependence of in-medium N-N cross section 
and insensitive to the symmetry potential. In the present paper, the calculations show 
that the conclusion drawn in Ref. [12] about Q zz remains true for the reversible kinetic 
reactions with heavy projectiles on light targets, but it is found that Q zz increases 
slightly as increasing the neutron-proton ratio of the colliding system. Therefore, Q zz 
can be a good probe for extracting information on in-medium N-N cross section in HIC 
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in the beam energies ranging from above Fermi energy to I50MeV/nucleon[ 121 , while 
the correlation between N im f and iV c may serve as a complementary one. 



2 Theoretical Model 

The dynamics of intermediate energy HIC described by QMDt 13-141 contains three 
ingredients: density dependent mean field, in-medium N-N cross section and Pauli 
blocking. To describe isospin effect appropriately, QMD should be modified prop- 
erly: the density dependent mean field should contain correct isospin terms including 
symmetry energy and Coulomb potential, the in-medium N-N cross section should be 
different for neutron-neutron ( proton-proton ) and neutron-proton collisions, and the 
Pauli blocking should be counted by distinguishing neutrons and protons. In addition, 
the initial condition of the ground state of two colliding nuclei should also contain 
isospin information. In the present calculation, the ground state of each colliding nu- 
cleus is prepared by using the initial code of IQMD, according to its density distribution 
obtained from the Skyrme-Hatree-Fock calculation with the parameter set of SKM*t 15 l. 
The interaction potential is as follows, 

U(p) = U Sky + V c (l - t z ) + U sym + V yuk + U MDI + U Pauli (I) 

jjSky j g density-dependent Skyrme potential, 

u Sk y = a (£~) + (2) 

Po Po 

V c is Coulomb potential. JJ Yuk is the Yukawa potential' 13 ^ , 

Uj uk = fsE*,- ^{e- r - /m [l - HVL/m - r«/2>/L)]- 

where $ is the error function. JJ MDI is the momentum dependent interaction ( MDI 

)[16] ) 

U MDI = Uln 2 [t^-n) 2 + ^\- (4) 

Po 

jjPauii ig the p auli potential 117 " 181 , 

uFm " = v " { ^r )3exf{Jli ^P l - M 2^ £ > i ™ < 5 > 

Polo z % A Po 



(3) 



1 for neutron-neutron or proton-proton 



for neutron-proton 
jjsym j g symmetry potential. In the present paper, three different kinds of JJ sym 
have been used' 1 ', 

Ur n = cF 1 (u)6T z (6) 



U* ym = cF 2 {u)5r z + l -cF 2 {u)5 2 
U s r i = cF 3 (u)5t z - \cF 3 (u)5 2 



(7) 
(8) 



with 



{1 for neutron 
— 1 for proton 

Here c is the strength of the symmetry potential, taking the value of 32MeV. F\(u) = u, 
F 2 (u) = u 2 and F 3 (u) = w 1/2 , u = 5 is the relative neutron excess 5 = e ^- £z = 22 ^ £jL - 
p and po, p n and p p are the nuclear density and its normal value, neutron density and 
and proton density, respectively The parameters of the interaction potentials are given 
in table I. 

Table 1. The parameters of the interaction potentials 
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P 


7 


h 


m 


u 


h 


v P 


Po 


<?o 


(MeV) 


(MeV) 




(MeV) 


(fm) 


(MeV) 


(MeV-' 2 ) 


(MeV) 


(MeV/c) 


(fm) 


-390.1 


320.3 


1.14 


7.5 


0.8 


1.57 


5 x 10" 4 


30 


400 


5.64 



It is worth mentioning that recent studies of collective flow in HIC at intermediate 
energies have indicated a reduction of the in-medium N-N cross sections' 19-21 '. An 
empirical expression of the in-medium N-N cross section' 21 ' is used: 

„med /i i „ P \—free i r\\ 

°NN = (1 + « — f^NN ( 9 ) 

Po 

with the parameter a ~ —0.2 which has been found to better reproduce the flow data 
[19-20] _ n eTe <7^. e is the experimental N-N cross section' 22 '. The neutron-proton cross 
section is about 3 times larger than the proton-proton or neutron-neutron cross section 
below 300 MeV. 
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Q zz will be used to describe the nuclear stopping: 

Q zz = E(2P,« 2 - PS? - Py{i) 2 )- (10) 

i 

In order to make the isospin effect on the multifragmentation process in HIC more 
prominent, comparable study is carried out for two pairs of reversible reaction systems. 
For each pair of comparable reaction systems the same mass of heavy projectiles and 
light targets, the same incident energy and the same impact parameter have been used 
as follows: ™Cd + 4 g Ar and ™Xe +20 Ca with the neutron-proton ratios 1.42 and 
1.16; lyZn+fl Ar and H Kr +2° Ca with the neutron-proton ratios 1.42 and 1.07. Hence 
the differences of N im f and Q zz for each pair of the comparable colliding systems are 
mainly due to the isospin effect on the process of multifragmentation and dissipation. 

By means of the modified isospin-dependent coalescence model' 23 ', we construct 
clusters within which the particle relative momentum is smaller than po= 300MeV/c 
and the relative distance is smaller than R = 3.5fm. To avoid the nonphysical clusters, 
the restructured aggregation model' 24 ' is used until no nonphysical cluster is produced. 

3 Results and Discussions 

N im j and N c are calculated event by event with the charge number of IMF from 
Z=3 to 13 for the heavy colliding systems, and Z=3 to 8 for the medium mass colliding 
systems. 

3.1 The isospin effect on the process of multifragmentation 
and dissipation in HIC 

Fig.l shows the time evolution of N im f for the central reactions 76 Zn + 40 Ar (line 
1) and 76 Kr+ 40 Ca (line 2) at E=80 MeV/nucleon (left panel), and 120 Cd+ 40 Ar (line 
1) and 120 Xe+ 40 Ca (line 2) at £=100 MeV/nucleon (right panel). From Fig.l one can 
see that N im f for the two neutron-poor colliding systems 120 Xe+ A0 Ca and 76 Kr+ 40 Ca 
are larger than those for the neutron-rich systems 120 Cd+ m Ar and 76 Zn + 40 Ar. This 
difference are mainly due to the isospin effect on the multifragmentation because other 
reaction conditions, except the ratios of neutron to proton, are the same. 
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As is well known that nuclear reaction mechanism and reaction yields are sensitive 
to both impact parameter and incident energy. Denoted by the same lines as in Fig.l, 
Fig. 2 depicts the multiplicity of intermediate mass fragments N im f as a function of 
impact parameter for the two pairs of colliding systems, H°Cd+fl Ar and ™Xe-\-^Ca, 
l^Zn+fl Ar and H Kr +20 Ca. The incident energy E is 100 MeV/nucleon for the heavy 
systems (left panel), and 80 Mev/nucleon for the medium systems (right panel). As 
in Fig.l, at small impact parameters, the lines labeled 2 for the neutron-poor colliding 
systems are always above the corresponding lines labeled 1 for the neutron-rich systems. 
The difference between them disappears gradually as increasing impact parameter. 

Fig.3 shows N im f as a function of beam energy from 15 MeV/nucleon to 200 
MeV/nucleon at 6=0 fm for the above two pairs of colliding systems with the line labels 
as in Fig.l. The relative locations between the two lines in each panel of the figure are 
also similar to those in Fig.l, namely, N im f for the two neutron-poor colliding systems 
120 Xe + 40 Ca and 76 Kr + 40 Ca are larger than those for the corresponding neutron- 
rich systems 120 Cd+ 40 Ar and 76 Zn+ /L0 Ar for the beam energy above 50MeV/nucleon. 
However, as the beam energy decreases to below 50MeV/nucleon the collision dynamics 
is governed by both the mean field and the nucleon-nucleon collisions. In this case, the 
difference between the two lines of the colliding systems in each pair vanishes gradually. 

From Figs.l, 2, and 3 we can see that the intermediate mass fragment multiplicity 
Nimf for the neutron-poor colliding system are larger than that of the corresponding 
neutron-rich system with the same mass projectile on the same mass target and with 
the same entrance channel conditions except that the ratios of neutron to proton of 
the colliding systems are different. 

The correlation between N im f and N c for the central collision of the above two 
pairs of colliding systems with the same beam energies as above is shown in Fig. 4. It 
is seen from the figure that the correlation between Ni m f and N c for the two pairs of 
colliding systems displays a clear isospin effect, i.e., the N im f-N c correlation for the 
neutron-poor systems is different from that of the neutron-rich systems. Considering 
the total yield of N im f (namely the integral area of each curve in the figure), one can 
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reach the same conclusion as from Fig.l, Fig.2, and Fig.3. This means that the Ni m f for 
neutron-poor systems (solid line), on average, is larger than that for the corresponding 
neutron-rich systems (dash line). The difference comes mainly from the isospin effect 
on multifragmentation in intermediate energy HIC. 

The mechanism of the above fragmentation process can be understood as follows. 
Experimentally the free neutron-proton cross section is about three times larger than 
the free neutron-neutron or proton-proton cross section below 300 MeV. The mean 
N-N cross section is defined as: 

N n p<J n p -\- (N nn ~\~ Npp^cTpp (3N n p -\- N nn -\- Npp) . 1N n p (■\-\\ 

N + N + N ~ N + N + N app ~ [ ~N~ )<T { ' 

1 ' l np~ ly pp~ ly nn ly np ~ 1 * nn ~ ly pp ly 

where N np , N nn and N pp are the collision numbers for neutron-proton, neutron-neutron, 
and proton-proton, respectively, and N = N np + N nn + N pp . a np , a nn , and a pp are 
the free N-N cross sections for neutron-proton, neutron-neutron, and proton-proton, 
respectively. In general, a nn = a pp = a. Because the total collision numbers for two 
colliding systems with the same mass projectile and the same mass target are the 
same, and the number of neutron-proton collisions for neutron-poor colliding system is 
larger than that of the neutron-rich colliding system, the mean total N-N cross section 
< a > of neutron-poor system is thus larger than that of neutron-rich system. Due to 
the effect of Pauli blocking, the effective collision numbers become smaller. But after 
considering the Pauli blocking, the above conclusion remains unchanged. Therefore, 
the neutron-poor system will have more effective N-N collisions and lead to stronger 
compression-expansion, resulting in a large number of multifragmentation Ni m f for the 
neutron-poor system in comparison to the neutron-rich system in the above 50MeV 
region where the N-N collisions are dominant. 

In Fig. 5 is plotted the time evolution of the quadrupole of single particle momentum 
distribution Q zz for the reactions 76 Zn+ A0 Ar ( dash line ) and 76 Kr + 40 Ca ( solid line 
) ( bottom panel), 120 Cd+ i0 Ar ( dash line ) and 120 Xe+ 40 Ca ( solid line) (top panel) 
at E=80, 100, 150MeV/nucleon and b = O.Ofm. In the figure, Q zz for the neutron-poor 
system is always smaller than that of the corresponding neutron-rich system. Smaller 
Q zz indicates larger dissipation of the initial longitudinal collective motion into the 



internal chaotic motion and the subsequent thermalization of the system. 

3.2 A good probe and a complementary probe to in-medium 
N-N cross section 

Since the reaction dynamics of HIC is mainly governed by both nuclear EOS and 
in-medium N-N cross section, to understand the collision dynamics in details, both 
ingredients should be studied carefully. As is well known that the effects of both ingre- 
dients are usually mixed in the dynamics and the main uncertainty of the information 
about the nuclear EOS extracted from HIC is due to our poor knowledge of the N-N 
cross section in medium. If one can find an experimental probe which can distinguish 
the contribution of EOS from that of in-medium N-N cross section, that will be very de- 
sirable. In this paper we have found that Q zz may serve as a good probe to in-medium 
N-N cross section. 

In Fig. 6 is given Q zz as a function of neutron-proton ratio for seven colliding sys- 
tems 76 Kr + 40 Ca, 120 Xe + 40 Ca, 6A Ni + 40 Ar, m Kr + 40 Ar, 76 Zn + 40 Ar, 85 Ge + 40 Ar, 
u Ni+ 47 Ar with the neutron-proton ratios 1.07, 1.16, 1.26, 1.33, 1.42, 1.5 and 1.56 at 
£'=100MeV/nucleon and 6=0.0 fm for four cases: (1) the symmetry potential JJl ym be- 
ing employed and in- medium N-N cross section cr^p being isospin-dependent, namely, 
Ui ym +a tso ; (2) JJi ym and N-N cross section a™^ being isospin-independent, denoted by 
jjsym + a n 0iso . ^ jjsym &nd ^ed being i sosp i n -dependent, denoted by Uf m + a iso ; (4) 

Ul ym and cr™^ being isospin-dependent, denoted by U^ ym + a tso . In fig. 6, lines labeled 
1, 2, 3 and 4 correspond to the above four cases. 

It is clear to see that Q zz depends strongly on the isospin dependence of in-medium 
N-N cross section and weakly on the symmetry potential ( namely, line 1 is located 
near lines 3 and 4, but far away from line 2 in Fig.6), though Q zz increases slightly 
with increasing neutron-proton ratio of asymmetry colliding system. In this case, Q zz 
is a good probe for extracting information on the isospin dependence of in-medium 
N-N cross section. We have discussed this in more details in Ref.[12], but it is not 
clear in Ref.[12] how Q zz depends on the neutron-proton ratio of colliding system. In 
addition, we shall report a complementary probe for extracting the information on 
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in-medium N-N cross section: the correlation between the multiplicity of intermediate 
mass fragments N im f and the total number of charged particles N c , based on the fact 
that it is also sensitive to the in-medium N-N cross section and insensitive to the 
symmetry potential in the chosen energy region. To study the contributions to N im f 
from different ingredients separately, we consider four cases as the same as in Fig.6. 

In Fig. 7 is plotted the time evolution of Ni m f for two different asymmetric colliding 
systems, 76 Zn+ A0 Ar at E = 80MeV/nucleon and b = O.Ofm (left panel), 120 Xe+ m Ca 
at E = lOOMeV/nucleon and b = O.Ofm (right panel). It is noted that lines 1 are 
always located near lines 3 and 4, but far away from lines 2 with increasing colliding 
time. This implies that the isospin effect of in-medium N-N cross section on N im f is 
more important than that of the symmetry potential in the energy region studied here. 

As is well known that nuclear reaction products and reaction mechanism sensitively 
depend on impact parameter. Fig. 8 shows the multiplicity of intermediate mass frag- 
ments Ni m f as a function of impact parameter for the four cases as illustrated in Fig.6 
for the reactions 76 Zn+ 40 Ar at E = 80MeV/nucleon ( left panel ) and l2Q Xe+ AQ Ca at 
E = lOOMeV/nucleon ( right panel ). As in Fig.6, lines 1 in the figures locate near lines 
3 and 4, but far away from lines 2 at small impact parameters. This indicates again 
that the isospin effect of in-medium N-N cross section on N im f is larger in comparison 
with that of symmetry potential. With increasing impact parmeter this isospin effect 
disappears gradually. 

The correlations between Ni m f and iV c for the four cases are plotted in Fig. 9 for 
the reactions 76 Zn + 40 Ar at E = 80MeV/nucleon ( left panel ) and 120 Xe + 40 Ca at 
E = lOOMeV/nucleon ( right panel ). The behaviors of the four kinds of lines in each 
panel are similar to those in Fig.6, namely, the integral area of the curve 1 is always 
close to those of the curves 3 and 4, but larger than that of the curve 2. From Fig. 7, 
Fig.8, and Fig. 9 one can see that the isospin effect of in-medium N-N cross section 
on the correlation between N^ m f and iV c is more important than that of symmetry 
potential. Here it should be stressed that we have to choose the incident energy for 
certain asymmetric colliding system carefully to get the above feature, namely, in this 
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energy region the N-N collisions should be dominant. In this case we may conclude 
that the correlation between the multiplicity of intermediate mass fragments N im f and 
the total number of charged particles N c can be used as a complementary probe to the 
isospin dependence of the in- medium N-N cross section in HIC. 

4 Summary and conclusions 

Starting from the simulation of the intermediate energy HIC by using IQMD, the 
calculated results have shown prominent isospin effects on the process of multifrag- 
mentation and dissipation, i.e., the intermediate mass fragment multiplicity N im f for a 
neutron-poor colliding system is always larger than that of the corresponding neutron- 
rich system, while the quadrupole of single particle momentum distribution Q zz for 
a neutron-poor system is smaller than that of the neutron-rich system with the same 
masses of projectile and target, and the same entrance channel conditions for small 
impact parameters. We also can see that Q zz increases slight with increasing ratio of 
neutron to proton in the colliding system. 

From the theoretical simulation, it is clear to see that Q zz depends strongly on the 
isospin dependence of in-medium N-N cross section and weakly on the symmetry po- 
tential in the beam energies ranging from about 50 MeV/nucleon to 150MeV/nucleon. 
And the correlation between N im f and N c has the same properties as Q zz in the chosen 
energy region. So we would suggest that Q zz may serve as a good probe and N im f 
a complementary probe for extracting the information on the isospin dependent N-N 
cross section in HIC. 
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Figure 1: The time evolution of N im f for the central collisions 120 Cd + 40 Ar (line 1) 
and 120 Xe + 40 Ca (line 2) at E = lOOMeV/nucleon (right panel), 76 Zn + 40 Ar (line 1) 
and 76 Kr + 40 Ca (line 2) at E = 80MeV/nucleon (left panel). 
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Figure 2: The multiplicity of intermediate mass fragments N im f as a function of 
impact parameter for the reactions 120 Cd + 40 Ar ( line 1 ) and 120 Xe + 40 Ca ( line 2 
) at E — lOOMeV/nucleon ( left panel ), and the reactions 76 Zn + 40 Ar ( line 1 ) and 
76 Kr + 40 Ca ( line 2 ) at E = 80MeV/nucleon ( right panel). 
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Figure 3: The N im f as a function of incident energies from 15MeV/nucleon to 
200MeV/nucleon at 6=0 fm for the reactions l2Q Cd+ m Ar ( line 1 ) and 120 Xe + 40 Ca 
( line 2 ) at E = lOOMeV/nucleon ( left panel ), and the reactions 76 Zn + 40 Ar ( line 
1 ) and 76 Kr + 40 Ca ( line 2 ) at E = 80MeV/nucleon ( right panel ). 




10 20 30 40 50 60 10 20 30 40 



N N 



Figure 4: The correlations between N im f and N c for the reactions 120 Cd+ 4 ° Ar ( dash 
line ) and 120 Xe + 40 Ca ( solid line ) at E = lOOMeV/nucleon ( left panel ), and the 
reactions 76 Zn+ m Ar ( dash line ) and 76 Kr+ m Ca ( solid line ) at E = 80MeV/nucleon 
( right panel ). 
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Figure 5: The time evolution of the quadrupole of single particle momentum distri- 
bution Q zz for the reactions 76 Zn + 40 Ar ( dash line ) and 76 Kr + 40 Ca ( solid line ) ( 
bottom panel ), and the reactions 120 Cd + 40 Ar ( dash line ) and 120 Xe + 40 Ca ( solid 
line )( top panel ) at £=80,100 and 150MeV/nucleon and b = O.Ofm. 




Figure 6: The quadrupole of single particle momentum distribution Q zz as a function 
of the neutron-proton ratio for seven colliding systems 76 Kr + 40 Ca, 120 'Xe + 40 Ca, 
M Ni + 40 Ar, m Kr + 40 Ar, 76 Zn + 40 Ar, 85 Ge + 40 Ar, and 7A Ni + 47 Ar at E = 
lOOMeV/nucleon and b = O.Ofm for the four cases ( see text ). 
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Figure 7: The time evolution of Ni m f 
80MeV/nucleon and b = O.Ofm ( left panel 
and b = O.Ofm ( right panel ). Lines 1, 2, I 
Fig.6. 




b(fm) 



for the reactions Zr + Ar at E = 
), 120 Xe + 40 Ca at E = lOOMeV/nucleon 
, and 4 correspond to the four cases as in 




b(fm) 



Figure 8: The multiplicity of intermediate mass fragments N im f as a function of 
impact parameter for the reactions 76 Zn + 40 Ar at E = 80MeV/nucleon ( left panel ) 
and 120 Xe+ 40 Ca at E = lOOMeV/nucleon ( right panel ) in the four cases as in Fig.6. 
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Figure 9: The correlations between N im f and N c for the reactions 76 Zn + 40 Ar at 
E = 80MeV/nucleon ( left panel ) and 120 Xe + 40 Ca at E = lOOMeV/nucleon (right 
panel) for the four cases as in Fig. 6. 
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